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Abstract 



We consider the production of the open charm and J ftp mesons in heavy ion 
collisions at BNL RHIC. We discuss several recently developed pictures for 
J/ip production and argue that a measurement at RHIC energies is crucial 
for disentangling these different descriptions. 
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I. INTRODUCTION 

The production of charmonium states J/ip and ip' have been measured in nucleus-nucleus 
(A+A) collisions at CERN SPS over the last 15 years by the NA38 and NA50 Collaborations. 
This experimental program was motivated by a suggestion of Matsui and Satz ]IJ to use the 
J /ip as a probe of the state of matter created in the early stage of the collision (see also Ref. 
for a recent review). 

In the original picture [Q], it is assumed that J/ip particles are created in primary nucleon- 
nucleon (N+N) collisions when the charm-anticharm quark pairs are produced. The fraction 
(a few percents) of initial charm pairs which formed a J/ip state is subsequently reduced in 
A+A collisions because of the J / ip suppression by a factor which depends upon the details 
associated with the impact parameter and energy of the collisions. 

The statistical approach, proposed in Ref. ||, assumes that J/ip mesons are created 
at hadronization according to the available hadronic phase-space similar to other (lighter) 
hadrons. We call this a hadron gas (HG) model. 

Recently a picture of the J/ip creation via coalescence (recombination) of charmed quarks 
was developed within several model approaches [ffl-H]. Similar to the HG model ||, the 
charmonium states are assumed to be formed at the hadronization stage. However, they 
are produced as a coalescence of earlier created c and c quarks. Therefore, the J/ip yield 
and the total charm yield are connected to each other and differ in general from the HG 
equilibrium values. 

A principal quantity to be studied in our paper is the ratio of the the J/ip multiplicity 
to the total number of cc pairs, N C c, created in A+A collisions: 

In the standard picture of J/ip suppression, the number of primary J/ip particles is propor- 
tional to the number of directly produced cc pairs with invariant mass below DD meson 
threshold M. The fraction of the subthreshold pairs in the total number of decreases 
when the collision energy yfs grows. The number of initially produced J/ip particles is 
further reduced because of their inelastic interactions with nucleons from the colliding nu- 
clei. This effect becomes stronger when the number of participants increases. In addition, 
J/ip particles can be destroyed by the secondary hadrons (co- movers) created in the A+A 
collision. The density of these hadrons increases with both N p and \/s. The decreasing of 
the ratio R due to the above mentioned effects is known as 'normal J/ip suppression 7 . The 
formation of quark-gluon plasma (QGP) at large values of y/s and N p is supposed to lead 
to a sudden and strong drop of the ratio (H). This is called 'anomalous J/ip suppression'. 
Therefore, there is an unambiguous consequence of the standard picture when the 

collision energy per nucleon pair \/s and/or the number of nucleon participants N p increase, 
the ratio (P decreases. 

In models where the J/ip mesons are produced by recombination of charm quarks at the 
hadronization stage, the picture can be much different. For the large number of cc pairs 
Acc >> 1 expected at the RHIC energies, the charmonium states can be formed from a c 
and c which were initially created in different cc pairs ||. Therefore, the multiplicity of J/ip 
due to recombination of charm quarks should be roughly proportional to: 
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m) ~ ^ > ( 2 ) 

^tot 

where N tot is the total multiplicity of produced hadrons. In this case the ratio ([I]) may be 
constant or even increase (we call this e J/ip enhancement') with N p and/or y/s depending 
on the underlying physics which generates cc pairs and charmonium states. 

We consider the canonical ensemble formulation || of the statistical coalescence model 
(SCM) [BJ to describe the c and c recombination into the charmonium states at the hadroniza- 
tion stage. A canonical rather than the more familiar grand canonical formalism is required 
as we must properly account for the small numbers of charm quarks which are present. 
When large numbers of charm quarks are produced in each collision, these two formalisms 
become equivalent, but the grand canonical formalism fails for small numbers. As in the 
standard picture, the charm is created at the very early stage of A+A reaction via hard 
parton collisions. The J/ip formation is however very different in the standard picture of 
the J ftp suppression |]J and in the SCM. This leads to very different relations between the 
open charm and charmonium yields. 

In the analysis below we present also two possibilities of charm chemical equilibration. 
First, we assume that both the open charm and charmonium yields correspond to a chem- 
ically equilibrated hadronic gas. This scenario is the most extreme in its assumptions, but 
provides a benchmark for the dispersion of various theoretical results. Another possibility 
is that the number of cc pairs, which is an input for the SCM, equals to the equilibrium 
value in the QGP just before its hadronization. These two possibilities are the extremes 
of scenarios with chemical equilibration of charm. We do not pretend that these extreme 
cases really represent what happens in heavy ion collisions. They are studied to complete 
our analysis. In fact, the hard parton collisions seem to be the dominant source of the cc 
pairs created at the RHIC energies. 

The issue we want to address in this paper is how to disentangle the different scenarios 
of charm and charmonium production. It will be shown that the different pictures have 
a much different dependence on both N p and y/s. Distinctive signatures of each of these 
pictures will be found. We shall argue that J / ip measurements at RHIC energies, especially 
in combination with a measurement of open charm, will go far towards resolving the various 
pictures. Of course, to really disentangle different scenarios, it is useful to be a little more 
quantitative. That is the central purpose of this paper. 

The outline of this paper is as follows. In Section 2 we present the equilibrium HG model 
for the open charm and charmonium production. The picture of the SCM is described in 
Section 3. The results and conclusions of the paper are presented respectively in Sections 4 
and 5. 



II. HADRON GAS MODEL 

The equilibrium HG model describes the hadron yields measured in A+A collisions in 
terms of three parameters: volume V, temperature T and baryonic chemical potential \xb- 
This model seems to work in the AGS-SPS energy region describing successfully not only 
the pion and (anti)nucleon yields but (multi) strange hadron abundances as well (see e.g. 
T0|). For the RHIC energies the temperature parameter T is expected to be similar to that 
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for the SPS energies: T = 175 ± 10 MeV. The baryonic chemical potential becomes small 
{Hb < T) and decreases with the collision energy. 

Does the HG model work for the charm and charmonium yields in A+A collisions? The 
answer on this question is rather uncertain at present. There are no open charm yield data 
measured in A+A collisions. This may be a reason why there were no attempts to calculate 
the open charm hadron yield within the HG model. 

A. Grand Canonical Ensemble Treatment of the Hadron Gas 

The HG model assumes the following formula for the hadron thermal multiplicities in 
the grand canonical ensemble: 



* = ^ L P'dP 



p 2 + mj- N 
exp | | ± 1 



(3) 



where V and T correspond to the volume^ and temperature of the HG system, mj, dj 
denote hadron masses and degeneracy factors. The particle chemical potential /ij in Eq.([3D 
is defined as 

/ij = bjfiB + Sj/j s + Cjfic , (4) 

where bj, Sj, Cj denote the baryonic number strangeness and charm of particle j. The baryonic 
chemical potential /i# regulates the baryonic density of the HG system whereas strange us 
and charm nc chemical potentials should be found from the requirement of zero value for 
the total strangeness and charm in the system (in our consideration we neglect small effects 
of a non-zero electrical chemical potential). 

Some of the particles which are observed after a nuclear collision result from the decay 
of resonances which have decoupled after the matter fell out of equilibrium. The total 
multiplicities Aj 0< of such observed particles in the HG model are therefore the sum of: 

Nf = Nj + Y,MR^J)N R , (5) 

R 

where By(R — > j) are the corresponding resonance decay branching ratios (both strong and 
electromagnetic decays are taken into account). The temperatures and chemical potentials 
used in these relationships are those at decoupling. The hadron yield ratios Nj ot /Nf ot in 
the grand canonical ensemble are therefore known functions of T and variables, and are 
independent of the volume parameter V. 

For the thermal multiplicities of both open charm and charmonium states the Bose and 
Fermi effects are negligible. This is because the charm quark mass is so large that the system 
is always dilute for charm quark densities. Therefore, Eq.(|3]) is simplified to: 

Nj = ^L__ T m) K 2 Uf j * d 3 V e^ T (^-J exp (-^fj , (6) 



L To avoid further complications we use ideal HG formulas and neglect excluded volume corrections. 
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where K2 is the modified Bessel function. The HG model gives the J/ip yield: 



Nj% = N m + Br(V')^' + Br( X i)iV Xl + Br( X2 )iV X2 , (7) 

where Nj/j,, N^, N Xl , N X2 are calculated according to Eq.(|6]) and Br(ip') ~ 0.54, Br(xi) ~ 
0.27, Br(%2) — 0.14 are the decay branching ratios of the excited charmonium states into 
J/V>. Note that T ~ 170H-180 MeV leads to the HG value of the thermal ratio of (if/) /(J /if;): 




0.04^0.05, (8) 



in agreement with data Jl I] in central (N p > 100) Pb+Pb collisions. This fact was first 
noticed in Ref. |]I2" |. At small N p as well as in p+p and p+A collisions the measured value 
of the (ip')/ (J/ip) ratio is several times larger than its statistical mechanics estimate 
Therefore, different versions of the statistical approach considered in this paper will be 
restricted to A+A collisions with N p > 100. We do not intend here to describe the open 
and hidden charm production in p+p, p+A and in very peripheral A+A collisions within 
the statistical models. 

Of course the temperature of decoupling of the hadron gas is quite large, and one should 
properly be worried about the consistency of such a description. It nevertheless does provide 
a fair description of available data on particle production. As such, it provides a standard 
against which interesting new physics might be calibrated. 



B. The Canonical Ensemble for a Hadron Gas 

The numbers of charm quarks produced in high energy nuclear collisions need not be 
large. Implicit in the grand canonical ensemble treatment is this large number assumption. 
The canonical ensemble treatment is valid for any small number of particles, and we now 
employ it to properly describe charm abundances. 

In the canonical ensemble formulation the thermal charmonium multiplicities are still 
given by Eq.(0) for a grand canonical ensemble as charmonium states have zero net charm. 
The multiplicities given by this equation (Q) for open charm hadrons are modified by an 
additional multiplicative 'canonical suppression' factor. This has been shown, for example, 
in Reference ||13||). The suppression factor is the same for all individual open charm states. 
Therefore, if Np is the total grand canonical ensemble multiplicity of all open charm and 
anticharm mesons and (anti)baryons, then the canonical ensemble (c.e.) value of the total 
open charm is equal to: 

M iV oJ 

where Io,h are the modified Bessel functions. To find No we use Eq.(^) for thermal 
multiplicities of the open charm hadrons in the grand canonical ensemble and take the 
summation over all known particles and resonances with open charm ||14|| . The factor 
Ii(No) I Iq{No) in this equation is due to the exact charm conservation. For N Q << 1 
one has I\ (No)/Iq(No) — Aq/2 and, therefore, the canonical ensemble total open charm 



5 



multiplicity Nq g - ~ Nq << N is strongly suppressed in comparison to the corresponding 
grand canonical ensemble result. At SPS energies the canonical ensemble suppression ef- 
fects, due to small numbers of produced charm, are important for the thermal open charm 
yield even in the most central Pb+Pb collisions. These suppression effects become crucial 
when the number of participants N p decreases. Note that for No << 1 the multiplicities 
of the open charm hadrons are proportional to V 2 in the canonical ensemble HG formu- 
lation (instead of V in the grand canonical ensemble). In the opposite case No » 1, 
Ii(No) I Iq(No) —>■ 1, therefore the grand canonical ensemble and canonical ensemble results 
coincide: Nq £ ' — > No- In a high energy collision of heavy nuclei, the average numbers of 
produced (anti)baryons, (anti)strange hadrons and charged particles are much larger than 
1. Due to this fact, the baryonic number, strangeness and electric charge of the HG system 
can be treated according to the grand canonical ensemble. This is not true in general for 
charm. 

The number of cc pairs in the HG, N^ G , is given by 

n " g = \ n °iB) + n »' (io) 

where Nh is the total HG multiplicity of particles with hidden charm. 

We can now calculate the J/ip multiplicity, since this involves only knowing the grand 
canonical ensemble results for cc bound states, and the total open charm multiplicity Nq € ^ 
since we now have the canonical ensemble results in hand @). We need the chemical freeze- 
out parameters V,T,fi B for A+A collisions. For the RHIC energies the temperature pa- 
rameter T is expected to be similar to that for the SPS energies: T = 175 ± 10 MeV. 
To fix the unknown volume parameter V and baryonic chemical potential /is we use the 
parametrization of the total pion multiplicity [|l5|j : 



N p ~ ° (v/i)V4 ' W 

where C m 1.46 GeV _1//2 and is the nucleon mass. Eq . (|TT|) is in an agreement with 
both the SPS data and the preliminary RHIC data in Au+Au collisions at y/s = 56 GeV 
and \fs = 130 GeV. The pion multiplicity ( |TT1) should be equated to the total HG pion 
multiplicity N^ ot calculated according to Eq.(0) which accounts for resonance decays into 
pionic final states. The HG parameters V and /is are found then as the solution of the 
following coupled equations: 

(tt) =N?(y,T,n B ) = Vu^T^b), (12) 
N p = Vn B (T,fi B ) , (13) 

where n B is the HG baryonic density. In these calculations we fix the temperature parameter 
T. 

The baryonic chemical potential for Au+Au collisions at the RHIC energies is small 
(fiB < T) and decreases with collision energy. Therefore, most of the thermal HG multiplic- 
ities become close to their limiting values at \i B —>■ 0. Consequently most of hadron ratios 
Nj 0t /Nt ot become independent of the collision energy. 
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If we assume that the charm quarks reach their chemically equilibrated values, then the 
chemical potential for charm vanishes. We can then compute the abundances of charmo- 
nium. The ratio of J/ip to negatively charged hadrons is only weakly dependent on the 
collision energy as a consequence of the considerations of the previous paragraph and equals 
approximately to: 



(h-) 



(0.6 1.5) 1(T 6 , (14) 



for T = 170 MeV and T = 185 MeV, respectively. The result (14) is close to the estimate 



H|l6j of this ratio, 1.0 10 6 , extracted from the NA50 data in Pb+Pb collisions at 158 



A GeV. 



C. Summary of Results 



In this section, we developed the grand canonical ensemble and canonical ensemble for- 
malism for describing abundances of particles in thermal equilibrium. Nucleon and pion 
abundances were used to constrain the parameters of these treatments: volume and baryon 
chemical potential. The temperature of the chemical freeze-out was assumed to be equal 
to its value at SPS energies. Formulae were derived for the abundances of open charm and 
hidden charm (cc states). In particular, the charm to negative charge particle multiplicity 
comes out correctly within such a naive treatment. 



III. STATISTICAL COALESCENCE MODEL 

One has serious reservation about the concept of charm chemical equilibration within HG 
approach. The HG system is expected to exist for a much shorter time than is needed for 
the chemical equilibration of charm. It is more reasonable to assume that the charm quark- 
antiquark pairs are created at the early stage of A+A reaction and the number of cc pairs 
remains approximately unchanged during the further evolution. After the hadronization, 
the system would correspond then to the HG picture with a nonequilibrium value of the 
total charm. 

We imagine two possible scenarios. In the first, we imagine the cc pairs reach their 
equilibrium values in the quark gluon plasma which precedes the hadronic gas. We call this 
the QGP + SCM model, where the SCM means statistical coalescence model. In the second, 
we imagine the cc pairs are produced by hard QCD processes early in the collision, and are 
unchanged by the presence of the hadronic gas or of the quark gluon plasma. This scenario 
we refer to as the pQCD + SCM model. 



A. Charm Quarks Not In Chemical Equilibrium 

The SCM |)]J| assumes that the charmonium states are formed at the hadronization 
stage. This is similar to the HG model [0]. However, the charmonium states are produced 
via a coalescence of earlier created cc quarks. The number of cc pairs, N C c, differs in general 
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from the result expected in the equilibrium HG. One needs then a new parameter 7 C || 
to adjust the thermal HG results to the required number of fixed by the early stage 
evolution. In the grand canonical ensemble, 7 C = e^/ T , where fi is the chemical potential 
for c + c pairs. This is analogous to the introduction of the strangeness suppression factor 
7 S |I7| in the HG model, when the total strangeness observed is smaller than its chemical 



equilibrium value. 

At high collision energies we find 7 C > 1 so that the open charm hadron yield is enhanced 
by a factor 7 C and charmonium yield by a factor 7^ in comparison with the equilibrium HG 
predictions. The canonical ensemble formulation of the SCM, as shown in Ref. 0, gives 

AT 1 AT ^{icNo) 2 AT / 1P -\ 

"* = r« N ° + 1- N " • (15) 

Note that the second term in the right-hand side of this equation is proportional to the 
number of hidden charm states, Nh, and is typically very small compared to the total 
number of charm quarks. It gives only a tiny correction to the first term. Most of the 
created cc pairs are transformed into the open charm hadrons. 

If N& » 1, Eq.flOjD is simplified to §: N c - C = 7c N Q /2 + ^N H . This happens for 
central Au+Au collisions at upper RHIC energy yfs = 200 GeV. For lower RHIC energy 
y/s = 56 GeV the value of iV cc - could be close to (or even smaller than) unity so that the 
canonical ensemble suppression effects for the open charm are still important. Note that for 
the non-central A+A collisions the total charm decreases essentially so that the canonical 
ensemble suppression effects become stronger, hence their consideration is necessary to study 
the N p dependence of the charmonium production even at the upper RHIC energy. 

The above equation for the number of cc pairs will be used to find the charm enhancement 
factor 7 C and calculate the J/ip multiplicity: 

m) = ii n i j% . (i6) 

where Nj°L is given by Eq.(^), the equation which accounts for decays of states with hidden 
charm into the J ftp. 



B. Estimating the Initial Charm Production 



The number of cc pairs, N C c, in the left-hand side of Eq.(|l5|) will be estimated in two 
different ways - in the pQCD approach, N c ^ = N^ CD , and in the equilibrium QGP just 
before its hadronization, N cZ = N^ GP , - and used then as the input for the SCM. 

The pQCD calculations for cc production cross section a(pp —>■ cc) in p+p collisions 
were first done in Ref. [0. In Ref. [fn| the cross section a(pp — > cc) was calculated in 
the next-to-leading order pQCD approximation and fit to the existing p+p data for charm 
production. We have parametrized the i/s-dependence of the charm production cross section 
at 10 < yfi < 200 GeV as: 




a(p P ^cc) = a It , (17) 
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with a ~ 3.392 fib, M ~ 2.984 GeV, a ~ 8.185 and (3 ~ 1.132. For our fit we used 
the result of Ref. jnj obtained with GRV HO set of proton structure functions |20| for 
charm quark mass and renormalization scale m c = fi = 1.3 GeV. The formula (|TT| ) agrees 
approximately with existing data and leads to the value of a(pp — > cc) ~ 0.35 mb at 
Vs = 200 GeV. 

The number of produced cc pairs in A+A collisions is proportional to the number of 
primary nucleon-nucleon collisions: 

N§ cr> = N A fi{Np) ^VV^f) , (18) 

where N^(N P ) is the number of primary nucleon-nucleon collisions, which depends on the 
number of participants and on the details of nuclear geometry, a™^ ~ 30 mb is the inelastic 
nucleon-nucleon cross sections. 

For the most central A+A collision (N p ~ 2A) the number of primary nucleon collisions 



can be estimated within the model of homogeneous spherical nucleus f21 |: 

*££{Np) * l^T (19) 

with Ra ~ 1.1 2 A 1 / 3 fm. Therefore N p dependence on the number of produced cc pairs can 
be calculated as:0 

N PQCD „ c a(pp ^ c _ } ^4/3 ^ (2Q) 

where C ~ 11 barn -1 . 



C. Initial Charm from a Chemically Equilibrated QGP 

We also make the SCM model calculations with N C c = N^- GP in the left-hand side of 
Eq. (|T5|) . The values of N® GP in central Au+Au collisions at different a/s are calculated as the 
equilibrium number of cc pairs in the QGP just before its hadronization. The hadronization 
temperature of the QGP is taken to be equal to the chemical freeze-out temperature in the 
HG, i.e. we assume that the chemical freeze-out in the HG takes place just after the QGP 
hadronization. The volume Vqgp of the QGP is found then from the requirement of an 
approximate entropy conservation during the hadronization transition: Sqgp = Shg (the 
nonzero but small baryonic chemical potentials in the HG and in the QGP do not affect 
the QGP volume estimate). The calculation of N® GP corresponds to the equilibrium ideal 
gas of quarks and gluons (m u « « 0, m s Pd 0.12 GeV, m c pa 1.3 GeV). In the grand 
canonical ensemble, one finds for the number of c and c: 



2 Strictly speaking, Eq.(|20|) does not exactly correspond to the real experimental situation. Au+Au 
collisions with different centrality rather than central collisions of nuclei with different atomic 
numbers will be studied at RHIC. Still, the formula (^) gives almost the same numerical result 
as the more realistic calculations 
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N c = N- c = 2S^E T nf e K 2 (^) , (21) 

where d c = 6 is the degeneracy factor for the c-quark. Taking into account the exact charm 
conservation we obtain the number of cc pairs in the canonical ensemble: 



,QCP . 1 ,„ , »„ h(N c + N z ) 



NT -= 5 (M^. (22) 



IV. COMPARISON OF VARIOUS MECHANISMS 

In this Section we present the results of the SCM and the HG model. Two different 
inputs for N c - C in the SCM (|5]) are studied: Njg CD © and N% GP ©. We call these 



model formulations as pQCD+SCM and QGP+SCM, respectively. The HG model, discussed 
in Section 2, can be considered as a special case of the SCM when 7 C = 1. First, we 
consider central Au+Au collisions (N p ~ 2A) at different collision energies y/s. Keeping the 
temperature parameter T fixed we solve Eqs.( [T2"| , |TB| ) to find the hadronization parameters 



V and hb- Eq-(Q) is used as the parametrization of the y/s dependence for the total pion 
multiplicity. With HG parameters T, V and /ig we calculate the thermal multiplicities 
of individual open and hidden charm states in the grand canonical ensemble according to 
Eq. (H) and then find the grand canonical ensemble total open charm Nq and the total hidden 
charm Nh multiplicities taking the summation over all known open and hidden charm states. 
Using the calculated values of N Q and N H we proceed further with three different model 
formulations along the scheme presented in Table 1. 



A. Number of cc Pairs as Function of Energy 

The dependence of the number of cc pairs, N C c, on the collision energy is shown in Fig. [I] 
for different model formulations. The value of N^ CD (p0[) increases strongly with collision 
energy and at the highest RHIC energy it becomes several times larger than the thermal 
equilibrium value iVjp in the HG model (pTOT) . At the SPS energies y/s m 20 GeV the 
situation is however different: iVjp > Nf!-^^ . The value of N® GP defined by Eq. (p2|) 
seems to give an upper limit of the charm production in A+A collisions. The inequality 
N^ GP » N^ G is mainly because of a small mass of charm quark, m c ~ 1.3 GeV, in 
comparison to the masses of charmed hadrons (D, D*, etc.). 



B. J/ijj Multiplicity versus Energy 



The J/ip multiplicity is shown in Fig. |2| and the ratio of J/ip to negative hadrons in 
Fig. [| The experimental estimate of the J/ip to h~ ratio extracted from the NA50 data in 
Pb+Pb collisions at 158 A GeV is also presented in Fig. |3|. It seems to be in agreement with 
simple HG results (i.e. 7 C ~ 1). In the pQCD+SCM model the enhancement of the total 
charm N C c by a factor of about 3 is needed to explain the observed J /if} multiplicity 0. This 
prediction of the statistical coalescence model can be tested in the near future (measurements 
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of the open charm are planned at CERN). Such a comparison will require to specify also 
more accurately the (J/ip) data. A possible source of the open charm enhancement in A+A 
collisions at the SPS with respect to the direct extrapolation fl20|) from p+p to A+A may 
be the broadening of the phase space available for the open charm production due to the 



presence of the quark-gluon medium ||23| . This effect, however, is expected to be small at 
RHIC energies. 



C. J/ip Multiplicity Ratio to Open Charm 

The dependence of the J/ip to N C c ratio on the collision energy is shown in Fig. [|. In the 
HG model and in QGP+SCM, the ratio of the J/ip to the number of cc pairs is constant at 
high collision energies. This is because both the J / ip multiplicity and the number of cc pairs 
are proportional to the system volume V. In contrast, the pQCD+SCM model predicts an 
essential increase of the ratio by a factor of about 3 from the lowest to the highest RHIC 
energy. 

Results presented in Figs. correspond to the most central Au+Au collisions (N p ~ 
2 A). The dependence of the ratio R on the number of participating nucleons N p is shown 
in Fig. |5]. (Note that statistical models are expected to be valid for N p > 100.) In the 
pQCD+SCM dependencies, one observes both the J/ij) suppression and the J/ip enhance- 
ment effects: the ratio (0) decreases with N p when the number of created cc pairs is small, 
N^ CD < 1, but it begins to increase with N p when N^ CD is essentially larger than unity. In 
the HG model and in QGP+SCM, the suppression is also observed at small N p corresponding 
to N^ CD < 1. In contrast to pQCD+SCM, the models with statistical charm production 
do not demonstrate J/ip enhancement: at large N p , which correspond to N^ CD » 1, the 
ratio R tends to a constant value. 



D. Centrality dependence of 7 C 

The dependence of charm enhancement factor 7 C on the number of participating nucleons 
N p is shown in Fig. |6|. In the case of pQCD+SCM model, the behavior is similar to that 
of ratio R: 7 C decreases with N p provided that the number of cc pairs is small N^ CD < 1 
and increases at large N^ CD . In contrast, the QGP+SCM model leads to a nearly constant 



value of 7c, which can be easily explaned compareing Eqs. ([15|) and (|22|) and neglecting the 
small term ^Nh- From the same reason the dependence of 7 C on ^fs is also negligeble. (Let 
us remind that in the case of HG model, 7 C = 1 by the definition). 



E. Discussion 



The results presented in Figs. [§-|6| dependences can be easily studied analytically in the 
limiting cases. For N C c >> 1 one finds from Eq.flT5|): 

R = ^ ~ 7 C - 2N C - C /N ~ ^ ~ (23) 
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where (ir) is the total pion multiplicity. Therefore, in the limiting case N C c » 1, the ratio of 
the J / ip to the open charm (or a probability to form the charmonium states) is proportional 
to the factor 7 C , which in its turn is proportional to the density of cc pairs. In the case of 
pQCD+SCM, this cc density increases both with collision energy y/s and with the number 
of nucleon participants N p : 

1) the increase (|20|) of iV^ D with ^/s is much stronger than the increase of the system 
volume V or, equivalently, the increase of the total pion multiplicity (pi]); 

2) the increase of N^ CD ~ N^ 3 is stronger than V ~ (n) which is approximately propor- 
tional to N p . 

In the case of thermal charm production (QGP+SCM and HG), the number of cc pairs is 
proportional to the system volume. Therefore the ratio R becomes constant when N C c » 1. 

The behavior of the (J /ip) /N& ratio is rather different when N c ^ « 1 and the canonical 
suppression factor in Eq . (|l5l) plays a crucial role. One finds 7^ ~ 4N C c/Nq and then j c ~ 

— 1/3 

N P for pQCD+SCM. The value of R, however, does not depend on 7 C : 

B= m.„L„±. (24) 

The ratio ([!]) decreases with both the collision energy (approximately like (^fs) 1 / 2 ) and 
the number of nucleon participants (approximately like N~ l ) and does not depend on the 
charm production mechanism. This behavior is similar to the standard picture of the J/ip 
suppression. 

The J/ip suppression in pQCD+SCM happens at the SPS energy. This energy is still 
"low" as N^ CD < 1 even in the most central Pb+Pb collisions. The J/ip suppression (i.e. 
the decrease of the J /ip/N v 9 CD ratio) with N p and/or s/s continues up to the RHIC energies 
for the peripheral collisions with small number of nucleon participants. This changes into 
the J/ip enhancement when both N p and \fs are "large". Therefore, if the cc production is 
described by the pQCD, the statistical coalescence model predicts the decrease of the ratio 
R ([!]) with y/s and/or N p (the J/ip suppression) when N^ CD << 1. However, this ratio 
increases with y/s and/or N p (the J/ip enhancement) when N^ CD » 1. These limiting 
behaviors are smoothly connected in the intermediate region 

Other models for J ftp production via coalescence of c and c at the hadronization 
stage ||,|7]] give similar qualitative result: the J/ip enhancement for central collisions 
at high energies. Quantitative predictions, however, are rather different. In Ref. 0] 
(J/ip) = (2.0 -T- 3.0) 10 _1 per central Au+Au event at s/s = 200 GeV, whereas our result 
in pQCD+SCM is (J ftp) — 1.0 10~ 2 (note that N C c estimates are very close: approximately 
10-=- 11 cc pairs per event). The predictions of the microscopical coalescence model of Ref. 
can be compared with our results at = 130 GeV. They are rather close for J / ip to nega- 
tive hadron ratio: (J/ip)/(h,-) = 1.6 10~ 6 and (J/ip)/(h~) ~ 2 10~ 6 in the pQCD+SCM. 
However, the total charm production is rather different: N cZ ~ 12 in Ref. [0] and N& ~ 5 in 
our pQCD estimation. 



V. CONCLUSION 

The statistical production of J/ip mesons in Au+Au collisions at RHIC energies has 
been studied. The obtained results are qualitatively different from the standard picture of 
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the J/ip suppression in A + A collisions |],|2|], if the average number of produced charmed 
quark- ant iquark pairs per event is large N C c » 1. The standard picture leads to monotonic 
decrease of the ratio R (the ratio of the J / ip multiplicity to the total number of charm quark- 
antiquark pairs) with both the collision energy and the number of participating nucleons. We 
predict another behavior. The reason for this difference is the completely distinct physical 
pictures for charmonia formation. The standard picture supposes the charmonia states to be 
created exclusively in primary nucleon-nucleon collisions at the early stage of the reaction 
so that all subsequent interactions can only destroy them. We study another possibility: 
creation of charmed and charmonium particles at the hadronization stage. In our opinion, 
the most probable scenario of the charmonium production in A+A collisions at the RHIC 
energies is the following (we call it pQCD+SCM): cc pairs are produced exclusively in 
hard parton collisions at the early stage of A+A reaction (pQCD), and these charmed 
quarks and antiquarks are redistributed in the open and hidden charm hadrons according 
to the statistical mechanics prescriptions at the hadronization stage (SCM). This model 
predicts a strong increase of the ratio R with the collision energy or the number of nucleon 
participants (J/ip enhancement), provided that the total number of produced cc pairs per 
event is essentially larger than 1. This condition is fulfilled at the highest RHIC energy 
\/i = 200 GeV already for N p = 100 and also at lower RHIC energies for the most central 
Au+Au collisions (N p ~ 2 A). 

To get a feeling of possible effects of the charm creation (annihilation) by secondary 
parton and/or hadron rescatterings we have considered two extreme models of charm chem- 
ical equilibration: in the quark gluon plasma (QGP+SCM) and in the hadron gas (HG). 
Both models predict almost constant ratio R at N C c » 1. There is however a quantitative 
difference: QGP+SCM predicts much larger value of R and much larger value of the total 
charm than pQCD+SCM and HG models. The predictions of the HG model for the open 
charm and charmonium production at the upper RHIC energies are below the results of the 
pQCD+SCM and QGP+SCM. 

Our predictions for the behavior of R at N C c » 1 {J/tp enhancement) in the 
pQCD+SCM are in a qualitative agreement with other coalescence models fH?], but quan- 
titatively the results are rather different. 

We conclude that Au+Au collision experiments at RHIC offer a possibility to disentangle 
various models of the open charm and charmonia production. Even those models that give 
very similar results at SPS energies become quantitatively and even qualitatively different at 
RHIC. Measurements of the J/ip productions combined with the open charm measurements 
are needed to give a definite answer about the dominant mechanisms of hidden and open 
charm production in ultrarelativistic collisions of heavy nuclei. 
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TABLES 
TABLE I. 





pQCD + SCM 


QGP + SCM 


HG 


N C c 


N^ CD , Eq. © 


N% GP , Eq. © 


iV c f G , Eq. ©) 


7c 


Eq. © , JV« = A^ QCD 


Eq.(|D , = A« GP 


7c = 1 




7c 2 i\% , Eq. (ID 


7c 2 i\% , Eq. (16) 
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FIGURES 




FIG. 1. The dependence of the average number of cc pairs per Au+Au central collision on the 
center-of-mass energy of the nucleon pair. The calculations are done in the framework of three 
different models: perturbative QCD (pQCD), equilibrium quark-gluon plasma (QGP), equilibrium 
hadron gas (HG). The QGP and HG results are shown for two different temperatures: T = 170 
MeV and T = 185 MeV. 
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FIG. 2. The dependence of the average number of J/^-particles per Au+Au central collision on 
the center-of-mass energy of the nucleon pair. The calculations are done in the framework of three 
different models: statistical coalescence model with the number of cc pairs given by perturbative 
QCD (pQCD+SCM), statistical coalescence model with the number of cc pairs given by equilibrium 
quark-gluon plasma (QGP+SCM) and equilibrium hadron gas model (HG). All the dependences 
are shown for two different temperatures: T = 170 MeV and T = 185 MeV. 
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FIG. 3. The dependence of the ratio of the J/ ^-multiplicity to the negative hadron multiplicity 
in Au+Au central collisions on the center-of-mass energy of the nucleon pair. The calculations are 
done in the framework of three different models: statistical coalescence model with the number of 
cc pairs given by perturbative QCD (pQCD+SCM), statistical coalescence model with the number 
of cc pairs given by equilibrium quark-gluon plasma (QGP+SCM) and equilibrium hadron gas 
model (HG). All the dependences are shown for two different temperatures: T = 170 MeV and 
T = 185 MeV. The filled circle shows the value estimated from the NA50 experimental data for 
Pb+Pb collisions at the SPS energy [2,14]. 
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FIG. 4. The dependence of the ratio of the J ftp- multiplicity to the multiplicity of cc pairs in 
Au+Au central collisions on the center-of-mass energy of the nucleon pair. The calculations are 
done in the framework of three different models: statistical coalescence model with the number of 
cc pairs given by perturbative QCD (pQCD+SCM), statistical coalescence model with the number 
of cc pairs given by equilibrium quark-gluon plasma (QGP+SCM) and equilibrium hadron gas 
model (HG). All the dependences are shown for two different temperatures: T = 170 MeV and 
T = 185 MeV. 
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FIG. 5. The dependence of the ratio of the J/Y>-rnultiplicity to the multiplicity of cc pairs 
on the number of nucleon participants N p for Au+Au collisions. The calculations are done in the 
framework of three different models: statistical coalescence model with the number of cc pairs given 
by perturbative QCD (pQCD+SCM), statistical coalescence model with the number of cc pairs 
given by equilibrium quark-gluon plasma (QGP+SCM) and equilibrium hadron gas model (HG). 
In pQCD+SCM case, the dependences are shown for three collision energies per the nucleon pair: 
RHIC-1 (y/s = 56 GeV), RHIC-2 {yfi = 130 GeV) and RHIC-3 (y/s = 200 GeV). For QGP+SCM 
and HG, the results at RHIC-1 and RHIC-3 are presented. The freeze-out temperature was fixed 
at T = 170 MeV. The models are assumed to be valid for iV p > 100, i.e. to the right from the 
dotted vertical line. 
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FIG. 6. The dependence of the charm enhancement factor 7 C on the number of nucleon par- 
ticipants N p for Au+Au collisions. The results are presented for three different models: statistical 
coalescence model with the number of cc pairs given by perturbative QCD (pQCD+SCM), sta- 
tistical coalescence model with the number of cc pairs given by equilibrium quark-gluon plasma 
(QGP+SCM) and equilibrium hadron gas model (HG). In pQCD+SCM case, the dependences 
are shown for three collision energies per the nucleon pair: RHIC-1 (^/s = 56 GeV), RHIC-2 
(^/s = 130 GeV) and RHIC-3 (^/s = 200 GeV). The freeze-out temperature was fixed at T = 170 
MeV. For QGP+SCM, the value of "y c is shown at two values of freeze-out temperature T = 170 
MeV and T = 185 MeV (the dependence on the collision energy is negligeble). For the equilibrium 
hadron gas model (HG) 7 C = 1 by the definition. The models are assumed to be valid for N p > 100, 
i.e. to the right from the dotted vertical line. 
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